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A phenomenological study of quark-hadron duality in electron and neutrino scattering on nuclei 
is performed. We compute the structure functions Fi and xFs in the resonance region within a 
framework that includes the Dortmund-group model for the production of the first four lowest-lying 
QO ■ baryonic resonances and a relativistic mean-field model for nuclei. We consider four-momentum 
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transfers between 0.2 and 2.5 GeV . The results indicate that nuclear effects play a different 



role in the resonance and DIS region. We find that global but not local duality works well. In 
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the studied range of four-momentum transfers, the integrated strength of the computed nuclear 
structure functions in the resonance region, is considerably lower than the DIS one. 



I. INTRODUCTION 



O ; n 

Nearly forty years ago, Bloom and Oilman found [1[ that in electron scattering on protons the inclusive structure 
—i, function F% in the resonance region oscillates around the DIS scaling curve and, after averaging, closely resembles 



it. This phenomenon is one of the ways quark-hadron duality reveals itself in physical processes. Generally quark- 
hadron duality establishes a relationship between the quark-gluon description of a certain phenomenon, which is 



theoretically justified in the DIS region, and the hadronic description, which is more convenient at medium and 
low energies. Understanding duality is also essential when establishing relationships between exclusive and inclusive 
processes. For a recent and detailed review of duality we refer the reader to Ref. [2|. 

So far, most theoretical studies of quark-hadron duality in lepton scattering were dealing with nucleon targets. The 



topic becomes of great practical interest when turning to nuclear targets and neutrino sources. The current precision 
measurements of the oscillation parameters require an efficient and accurate description of the neutrino-nucleus cross 
sections. Of particular interest is the resonance region and the possibility of linking it with the DIS region. A hadronic 
description of the neutrino-nucleus cross sections at low Q 2 requires the vector and axial transition form factors for 
each resonance. For the majority of the resonances, these transition form factors are not well constrained. Provided 
that one can establish that quark-hadron duality holds with a reasonable accuracy, one could think of using the DIS 
results for estimating the neutrino-nucleus cross sections in the resonance region. In that respect it is worth mentioning 
that in nuclei, the Fermi motion of the nucleons smears the observables, so that the averaging in the resonance region 
required for duality, proceeds to a certain extent automatically. The issue whether quark-hadron duality holds 
with sufficient accuracy in lepton-nucleus scattering, requires further theoretical and experimental investigation. The 
present paper addresses this issue from the theoretical point of view. 
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Recent electron scattering measurements at Jefferson Laboratory (JLab) have confirmed the validity of Bloom- 
Gilman duality for the proton, deuterium 3] and iron (jj structure functions. Further experimental efforts are required 
for neutrino scattering. Among the upcoming neutrino experiments, Miner^aj^, 3, Q] and SciBooNEflflQ aim at 
measurements with carbon, iron and lead nuclei as targets. From the theoretical side, recent investigations of the 
phenomenon of duality for electron and neutrino scattering on nucleons include the works reported in Refs. [ill. H. . 
These studies differ in the way they treat the resonant contributions and the way they parameterize the DIS structure 
functions. This paper extends the study of Ref. 13| about the duality phenomenon in the nucleon to nuclei. 

For a free nucleon target, the structure functions generally depend on the transferred energy v = E — E' and 
four-momentum Q 2 = —q^q^. At low Q 2 the ^-distributions reveal several peaks, which correspond to various baryon 
resonances. We briefly sketch our theoretical approach to resonance production in nuclei in Section [II Al The nuclear 
structure functions are defined in Section III Bl At high Q 2 the structure functions exhibit scaling behavior, which 
is discussed in Section III CI Comparing the structure functions in these two regions allows one to check the basic 
features of duality and compare its validity for different targets and incoming leptons. Our results are presented in 
Section [TTTI for electrons and Section ITVl for neutrinos. Conclusions are given in Section fVl 



II. FORMALISM 

We consider inclusive charged-current (CC) neutrino scattering from nuclei and its electromagnetic counterpart 

v^w) + a -> r(fc^) + x , ro^ + A^ro'^ + x , (i) 

where I is the lepton flavor, A represents a nucleus with mass number A, and k^ = (E, k) and k /fl = (E', k') are the 
four-momenta of the incoming and outgoing lepton respectively. We work in the laboratory frame of reference. The 
coordinate system is chosen such that the z-axis lies along the direction of the virtual photon, so that the transferred 
momentum is given by q^ L = fc' 1 — fc'' 1 = (v, 0, 0, q z ). The lepton scattering proceeds in the xz-plane. In this section, 
we investigate the structure functions Fq, , 2xF\ and xF% , the latter being nonzero for neutrino reactions only. To this 
end, CP-violation effects are neglected for the case of electron scattering. 



A. Resonance production on a nucleus 

For lepton-nucleus scattering we describe the struck nucleus as a collection of bound nucleons. Assuming an 
independent-particle shell model, each nucleon occupies a nuclear shell a with a characteristic binding energy e a 
and is described by the bound-state spinor u a . In the impulse approximation, an impinging lepton interacts with a 
single bound nucleon. Hence, the nuclear cross section can be expressed as an incoherent sum over all nucleons of 
onc-nucleon cross sections weighted with the corresponding nucleon momentum distributions n a . For example, for a 
carbon nucleus, one has 
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TABLE I: Binding energies (MeV) for carbon and iron nuclei 
proton neutron 



12 C: 






Is 1 / 2 


47.76 


51.17 


lp"' 2 


16.76 


19.87 


56 Fe : 






ls 1/2 


57.19 


63.66 


lp 


A^K 1 1 




Ip 1 ! 2 


39.32 


46.00 


Id 5 ' 2 


27.64 


34.84 


2s 1 ' 2 


17.77 


24.41 


Id 3 / 2 


16.55 


23.01 


lf /2 


12.11 


19.17 


2p 3 ' 2 




5.99 



This allows us to employ the one-body lepton-nucleon vertex that can be well constrained in experiments with a 
proton and deuteron target. The four-momentum of the bound nucleon can be written as — (m^ — e a ,p). Both 
the bound-state spinor u a (p) and the corresponding binding energies are computed in the Hartree approximation 
to the a — to Walecka-Serot model H] ■ Binding energies for carbon and iron are summarized in Table |TJ For 
each shell, the nucleon momentum distribution n Q (|p|) is constructed from the bound-state spinors, the normalization 
convention being 

d 3 p n a (\p\) = 1 . 

These n- Q (|p|) are shown in Fig.[Tl for the case of a carbon nucleus. Clearly, for a specific shell, the proton and neutron 
distributions are almost identical. 
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FIG. 1: Momentum distributions for proton and neutron shells in carbon. 



After the interaction takes place inside the nucleus, the reaction products can escape the nucleus without interactions 
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or they can undergo elastic and/or inelastic rescatterings with the other nucleons. Thus, the reaction strength is 
redistributed between different channels. All these processes are called the final state interactions (FSI). The effect 



of FSI can be large for a specific exclusive process, for example for quasi-elastic nucleon knockout [lq |. where the 
cross section can be suppressed by a factor of 2. In one-pion production, the outgoing pion can be absorbed in the 
nucleus and thus mimic a quasi-elastic event. For a duality study, however, it suffices to consider inclusive reactions. 
Consequently, since the outgoing hadrons and the residual nucleus are not detected, we can make the assumption, 
following Ref. [17], that FSI can be disregarded. 



Recently, duality in lepton-nucleon scattering was investigated theoretically within the Sato-Lee [18], Rein Sehgal 



19 1 and Dortmund-group 20( models for resonance production. In this paper, we follow the approach used in [201 ] and 
extend it to calculate the nuclear structure functions. In particular, in the resonance region we take into account the 
first four low-mass baryon resonances P33(1232), Pn(1440), Z?i3(1520), Sii(1535) and describe the vertices of their 
leptoproduction within a phenomenological form-factor approach. The nucleon structure functions Wj are defined by 
the standard expansion of the hadronic tensor 

W MV = -g^Wi + -^Wz - le^xv—^-Ws + -^Wi + — § Ws • ( 2 ) 

Each Wj depends on two independent kinematic variables, for example Q 2 and v, which are determined exclusively 
by the lepton kinematics. Another set of variables, namely Q 2 and W, is also possible, since the invariant mass W, 
defined as W 2 = (p + q) 2 , for a free target nucleon can be uniquely related to Q 2 and v. W 2 = m 2 N + 2m^v — Q 2 . 
The analytical expressions for the one-nucleon structure functions F± — m^Wi, F2 — VW2, F3 = 1SW3 in terms of 
form factors for a free nucleon as well as the form factors themselves are given in 20]. The Fermi motion of the 
bound nucleon modifies the expression for the scalar product (q ■ p), so that the invariant mass W 2 — (p + q) 2 will 
now depend on the nucleon momentum and binding energy. The variables Q 2 and v, being determined by lepton 
kinematics only, remain unaffected. Strictly speaking, the expansion in Eq. ([2]) is only valid for a free (on-mass shell) 
target nucleon. For a bound nucleon, all inclusive observables depend not only on v and Q 2 , but also on an additional 
independent kinematical variable, which can be chosen to be p^p^ = p 2 ■ Here, we make the assumption (see 21] 



for a detailed discussion) that expression © can still be used to define the bound-nuclcon structure functions, and 
recalculate them keeping the kinematical variable p 2 as an independent one. The results are given below for the 



W2 and W3 structure functions. Equating p 2 = mfj, the free-nucleon results of 20( are easily reproduced. For the 
spin-3/2 resonances (P33(1232) and Di3(1520) in our case) one has 

Wi(Q 2 ,u,p 2 ) = -^^(Q 2 ,^,p 2 )i?(^,Af fl ), (3) 
3mjv 

where R(W, Mr) is the finite representation of the 5— function <5(VF 2 — Mi), which gives the relativistic Breit-Wigner 
distribution: 

R{W > Mr) = tt (W 2 - M 2 ) 2 + M 2 R T 2 R ' 
and the Vi are given below. The upper and lower signs are for the positive (P33(1232)) and negative (Di3(1520)) 
parity resonances, respectively. 
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For spin- 1/2 resonances we have 



[2M% T 2m N M R + Q 2 - q ■ p] + 2 
2M% ± 2m N M R + Q 2 - q ■ p) 



c v b c 5 



(Q 2 -q-p) 
CYCt 



V r<A w 5 w 4 (rfl 



(Q 2 -q-p) 



N 



q-p 



(5) 



Wi(Q 2 ,v,p 2 ) = Vi(Q 2 ,u,p 2 )R(W,M R ) 
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and fi = tojv + Mr. The upper and lower signs again correspond to positive (Pn(1440)) and negative (511(1535)) 
parity resonances, respectively. In the case of electroproduction, all axial form factors should be put equal to zero 
and the weak vector form factors should be replaced by the electromagnetic ones for proton or neutron, depending 
on the target nucleon. To make the article self-contained, we present the transition form factors for each resonance. 
Electromagnetic and weak vector form factors were determined in [20I bv fitting the electroproduction data on helicity 
amplitudes in the region Q 2 < 3 GcV 2 . Recently, it was shown [22| that in order to satisfy the asymptotics for 
helicity amplitudes at Q 2 — *■ 00, as prescribed by perturbative QCD, the vector form factors should also exhibit a 
certain asymptotic Q 2 behavior. Therefore, we refitted the form factors according to this prescription. In the region 



Q 2 < 4 GeV 2 , however, the difference between our new fit and the one performed in 20 1 falls within the accuracy 
of the experimentally extracted helicity amplitudes. To be on the safe side for higher Q 2 values, further attempts 
to improve the fits of the form factors (for example, in accordance to upcoming data on helicity amplitudes) will be 
done within the framework of the arguments presented in [22j. The axial form factors are the ones used in 13 1 for 
the "fast" fall-off case. Thus, we use the following form factors 
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Here, D v = (1 + Q 2 /M 2 ) 2 with My = 0.84 GeV and D A = (1 + Q 2 /Af^) 2 with M A = 1.05 GeV. The weak form 
factors presented here are determined for the excitation of the R + resonance state, i.e. for neutrino scattering on a 
neutron. For the excitation of the double charged states, which is possible for isospin-3/2 resonances in neutrino- 
proton scattering, the isospin relation gives an additional factor y3 for each form factor. 



For the resonance widths we use the so called running widths Tr(W), as they were presented in Ref. 2fJ: 

2s R 



r R (w) = r° R 



pAw) 



(12) 



where sr is the spin of the resonance, on-shell widths are = 0.12 GeV, Fp 1440 = 0.350 GeV, r^, 1520 = 0.125 GeV, 
rsi535 = 0-150 GeV, and 

Pn{W) = ^\J( W2 - m2 N - m l) 2 ~ ^m 2 N ml . 



Definition of the nuclear structure functions 



For nuclear targets, the nuclear structure functions Wj can be defined in the standard manner by means of the 
expansion of the nuclear hadronic tensor 



W± = - gixv W? ' P » P »^ A 



[lis 



A A 

l£ ^ X °2Mi W * + Mf W4 + M\ Ws ' (13) 

where p A ' = (Ma, 0) is the four-momentum of the target nucleus with mass Ma in the laboratory frame of reference. 

In the impulse approximation we are dealing with the bound nucleon as a target, so we must relate the one-bound- 
nucleon structure functions introduced in the previous section to nuclear ones. We follow the prescription of Ref. 23 1 
and express the nuclear structure functions in terms of the nucleon ones in terms of a convolution of the type 
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(14) 
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where a extends over single-particle shells in the target nucleus and 2j a + 1 specifies their occupancies. 

It is worth stressing that in the original paper [23| as well as in 2l| an additional phase-space correction factor 
Ep/rriN is introduced in the expression (|14p to preserve the space volume under Lorentz transformation. Since we 
construct a momentum distribution from wave functions normalized as u^ a u a = 1 for each shell a, our correction 
factor must be equal to 1. 

Substituting ([2|) and (fT3| in (fT4| . one arrives at 
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This prescription guarantees, that as Q 2 tends to zero, the longitudinal structure function Wl also tends to zero as 
expected for the real photon: 

,2 



lim 



rW 2 A (Q>)-Wf(QV) 



. 



In neutrino experiments one can also measure the W3 structure function, for which our definition gives: 

Ma p°q z — vp z 
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Note that W3 depends on the nucleus mass Ma- Realizing that the Bjorken variable for a nucleus (xa = Q 2 /2Mav) 
differs from the one for a nucleon (x = Q 2 /2mjqv), the function that is independent of Ma is xaF%: 



x A F£ = E / d3 P ( 2 ^« + 1K(p)xF 3 (Q 2 , u,p) 



1 p°q z 
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m N q z 

According to the definition (fT"5"|) . F^ — an d xaF^ — xaMa^V\ are also independent on Ma- 

Within the adopted approach there is no unambiguous recipe for deciding whether one should keep mjv in the 
denominators of (|15[) or replace it with some effective mass, that corrects for the binding energy. For the numerical 
calculations presented here, we have opted to use the expression (|15p and interpret the m^r as the free nucleon mass. 

The integration over d 3 p = \p\ 2 d\p\ dcosj p dip p in Eq. (|15[) is performed in the following way. Integration over the 
azimuthal angle dip p gives 2-7T, since no structure function depends on it. The phase space in the plane determined 
by the absolute momentum value \p\ and polar angle j p is restricted by the condition W 2 > W^ in . For one-pion 
production one has that W m i n = + tn v . For a a bound nucleon this condition translates into 
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When performing the d 3 p integrations, the above condition determines the boundaries of the absolute bound-nucleon 
momentum \p\ for a given cos7 P , Q 2 and v: 



101= 
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The sign of the quantity W^ in + Q 2 —Pg — 2p§v discriminates between two classes of kinematic conditions. In what 
follows we provide a discussion of the values of p m in and p ma x in the phase-space integration JJ max dp for a positive 
and negative sign of W^ in + Q 2 — p\ — 2p§v. For 
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Q 2 -p 2 -2p v<0 , 



(21) 
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the calculated according to (|20[) is negative, so one should take p m in(Q 2 , v, cos7 P ) = 0. This means that the 
phase space (|21l) is accessible for a nucleon with arbitrarily small three-momentum, including \p\ = 0, as is the case 
for a free nucleon. When the condition (f2~Tj) is fullfilled, p m ax(Q 2 , v, cos7 P ) = |p| + (Q 2 , ^, cos7 P ) for all polar angles 
1p- 

Increasing the phase space for the bound nucleon does not necessarily imply that the cross section grows, because 
each point in the phase space gets weighted with a momentum distribution of the type shown in Fig. [TJ Cross sections 
and structure functions for high \p\ are strongly suppressed and the major contributions stem from the momenta 
inside the Fermi sphere. 

For 

W 2 mm +Q 2 -pl-2p»v>Q, (22) 

the \p\± are only defined for backward moving target nucleons. The restrictions on cos7 p for given Q 2 and v come 
from the condition 

(Q 2 + v 2 ) cos 2 7p + P 2 + 2 Po v -Q 2 ~ W 2 nm > , (23) 

which gives 

, , m „ (n i ,a . l w L n + Q 2 -pI- 2 p^ f9A s 

- 1 < cos7 p (y ,v) < -W q 2 + v2 • (^4) 

Since the minimal value of the three- momentum is positive in this case, the accessibility to this (Q 2 , v) region 
crucially depends on a nucleon already moving, which is only possible for a bound nucleon. This region of phase space 
grows in importance with increasing Q 2 . 

For Q 2 = 0.1 GeV 2 and different v, the typical phase spaces available for a Is 1 / 2 proton in carbon are shown in 
Fig. [2 We use polar coordinates for the variables \p\ and 7 P . The left (right) panel corresponds with the condition 
(E2) ((EH))- F° r each Q 2 and v, thick lines represent \p\ + and thin lines \p\~- The points where the \p\ + and 
lines coincide correspond to the upper boundary on cos7 P , as calculated in Eq. (|24p . Remark that the available phase 
space in (|p|,cos7 p ) is contained within a circle. At v — 0.4 GeV the left part of the circle is not shown because it 
corresponds to values of \p\ larger than 1 GeV. The momentum distribution of nucleons in nuclei will reduce those 
contributions to negligible proportions. 

The phase space collapses to one point cos7 P = — 1, \p\ = \/Q 2 + v^in f° r v min — ~Po + W m i n . Remark that for a 
bound nucleon the minimal value of v does not depend on Q 2 . Physically this means that for any Q 2 there is a bound 
nucleon moving backward fast enough to fulfill the requirement (jp + q) 2 > W 2 lin . Thus, contrary to the free nucleon 
case, for the off-shell nucleon the pion production threshold is defined in terms of v rather than invariant mass and, 
strictly speaking, is independent of Q 2 . 

At high Q 2 , however, using v m i n = —po + W m i n is not convenient for calculations, because all observables are 
strongly suppressed for large p inspite of the fact that the phase space is available. In our numerical calculations 
we have not considered nucleon momenta beyond three times the Fermi momentum. We stress that the phase space 
boundaries derived here depend on our assumption about the form of the four-momentum for the bound nucleon, 
which was taken as p^ = (mjv — & a ,P)- 
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FIG. 2: Sketch of the available (|p| , 7 P ) phase space in polar coordinates. We consider a ls 1//2 proton in carbon for Q 2 = 
0.1 GeV 2 and different v. The |p| is expressed in GeV. The left (right) panel corresponds with kinematics conditions obeying 
the condition of Eq. (J22j| (of Eq. (pTJ)). 



C. DIS region and scaling variable 



In the kinematical regime of high Q 2 and v, the so-called Bjorken limit, the structure functions depend only on 
the Bjorken variable x = Q 2 j2rajqv when one neglects higher-twist effects. This phenomenon of no observed Q 2 
dependence for a fixed x value is called Bjorken scaling. At these energies, the lepton scattering on nucleons and 
nuclei is dominated by deep inelastic scattering with a multiple-particle hadronic final state. Deep inelastic scattering 
on nuclei was intensively studied experimentally since the sixties. This experimental information will be used as DIS 
input for our investigation. 

For electron-carbon scattering, F2 was measured by the BCDMS Collaboration Q Q for 30 GeV 2 < Q 2 < 
200 GeV 2 . We choose several sets of data at different Q 2 — Q 2 DIS '- 30, 45 and 50 GeV 2 . As expected from Bjorken 
scaling, for most of the x region the data coincide with an accuracy better than 5%. For iron, the neutrino scattering 
results are available from the CCFR 26] and NuTeV [23| collaborations. 

Scaling structure functions are conventionally plotted against the Bjorken variable x. Violation of Bjorken scaling 
comes from target-mass corrections and higher-twist effects. In the scaling region, the Nachtmann variable £ — 



2x/ (1 + ^/l + Am 2 N x 2 / Q 2 ) was shown [28, 29] to be a better alternative, because it implicitly includes the kinematical 
part of the target-mass correction, which can be important at large x and low Q 2 . Expanding the inverse of this 
variable in a power series of 1/Q 2 , we recover the variable 1/^ « uj' — (^m^v + m 2 N )/Q 2 , used by Bloom and Gilman 
in their pioneering work on duality. For large v, one has lu' w 1/x. 



III. DUALITY IN ELECTROPRODUCTION 

In the case of an isoscalar target nucleon, and for Q 2 > 0.5 GeV 2 , it was shown [13j that Bloom-Gilman duality 
holds at the level of 20%. Here, we compute the nuclear structure functions and xaF\ along the model outlined 
in Section II. The results of our description of the resonance region in terms of hadronic degrees of freedom are then 
compared to DIS data. 
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In Figure [3l the bound-nucleon structure functions F2 and 2xF\ for a proton in the Is 1 / 2 and lp 3 / 2 carbon shells 
are contrasted with the structure functions for a free proton. They are plotted versus the Nachtmann variable £ for 
Q 2 = 0.2, 0.85, and 2.4 GeV 2 , with the largest Q 2 curves covering the largest £ values. Similar to the free-nucleon 
case, for a fixed Q 2 , the peak at larger £ corresponds to the A resonance and the peak at smaller £ corresponds to the 
second resonance region. One can easily notice the effect of smearing: the two resonance regions are distinguishable 
only at low Q 2 . Fermi smearing proceeds differently for different shells, which in turn introduces an additional 
averaging when summing over shells. One can also observe that the bound-nucleon curves extend to higher £ values 
than the free nucleon ones. This additional contribution comes from the phase space at low v values, which is shown 
in the left panel of Fig.[2]and discussed in Section Hi Bl At high £ the Fi and 2xF\ for the Is 1 / 2 shell are significantly 
lower than for the lp 3 / 2 shell. At high £ the phase space extends to relatively large bound-nucleon momenta. For 
those momenta the momentum distribution for a shell close to the Fermi surface (like lp 3 / 2 in carbon) is larger than 
for a deep-lying shell (like Is 1 / 2 in carbon). 
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FIG. 3: (color online) Structure functions F2 (left) and 2aj.Fi (right) for a free proton (solid curve), for Is 1 / 2 (dashed curve) 
and lp 3//2 (dash-dotted curve) protons in 12 C. The three sets of curves correspond to Q 2 = 0.2, 0.85, and 2.4 GeV 2 

Fig. [¥] shows the carbon structure function per nucleon F% c / A in the resonance region for several Q 2 values, 
from 0.45 to 3.3 GeV 2 . When investigating duality for a free nucleon, we took the average over free proton and 
neutron targets, thus considering the isoscalar structure function. Since the carbon nucleus contains an equal number 
of protons and neutrons, averaging over isospin is performed automatically. At Q 2 = 0.45 GeV 2 , the A peak is 
pronounced and can still be distinguished from the second resonance peak, which is also visible. At higher Q 2 one 
cannot distinguish the resonance structure anymore and the first and second resonance region merge into one broad 
peak. 



In the left panel of Fig. 0] the resonance structure functions are compared with data obtained by the BCDMS 

t2\ 



Collaboration 



24 



25j in muon-carbon scattering in the DIS region (Q 2 - 30 - 50 GeV^). They are shown as 



experimental points connected by smooth curves. For different Q 2 values, the curves agree within 5% in most of the £ 
region, as expected from Bjorken scaling. One observes that, as Q 2 increases, the resonance peaks decrease in height 
and slide along the DIS curve. This means that global duality holds for electron scattering on nuclei. To characterize 
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4: (color online) Duality for the F| c structure function. (Left) Resonance curves Ff c /12 as a function of £, for 



Q 2 = 0.45,0.85,1.4,2.4 and 3.3 GeV 2 (indicated on the spectra), compared with the experimental data [24J, [3] in the DIS 



region at Q DI s = 30, 45 and 50 GeV . (Right) Ratio I2 defined in Eq. (|25|) for the free nucleon (dash-dotted line), and c. We 
consider the under limits determined by W = 1.1 GeV (solid line) and by the threshold value (dotted line). 



local duality, we consider the ratio of the integrals of the resonance (res) and DIS structure functions 

C?max ,j t ^(res) . 



it ^<->\ 



(DIS), 



(25) 



>(f O 2 ) 

where Ti denotes or xaF^ (used later for neutrino scattering). The value Qjjjs is taken as the actual Q 2 value 
for a given experimental data set. For electron-carbon scattering we choose the data set at Q 2 DI s = 50 GeV 2 , 
because it covers most of the £ region. For a proton target [3[ , the integration limits for £ are conventionally chosen 
equal for both integrals and are defined in such a way as to cover the first and second resonance regions for each Q 2 . 
For a free nucleon, this requirement is written as 13 1 



£mi„ = £0^=1-6 GeV, Q 2 ), 



i{W = 1.1 GeV, Q 2 



(26) 



where the invariant mass for a free nucleon can be expressed in terms of v and Q 2 as W 2 — (p+q) 2 = m 2 y + 2mNis—Q 2 . 
The upper value W = 1.6 GeV is chosen in such a way as to cover the mass range of the four resonances taken into 
account, the heaviest one with the mass Mr — 1.535 GeV. The lower value W — 1.1 GeV is chosen close to the 
pion-production threshold Wthr = 1-08 GeV. In a nuclear target, the invariant mass of the struck nucleon depends on 
the initial momentum of the target nucleon. On the other hand, the structure functions, as well as other observables, 
are defined as integrals over the initial nucleon momentum. This prevents one from using W in defining the integration 
limits. One needs an alternative variable, which can be easily determined from the lepton kinematics. 

Experimentally one often (see, for example, "defines" the effective variable W by the relation W 2 = m 2 N + 

Ithnis — Q 2 . Notice that W is only an invariant for p = 0. However, it gives a reasonable feeling of the invariant mass 
region involved in the problem. In particular, the resonance curves presented in all figures are plotted in the region 
from the pion-production threshold up to W = 2 GeV. As was illustrated in Fig. [2l bound backward-moving nucleons 
allow lower v values beyond the free-nucleon limits. Thus, as discussed at the end of Secton III B[ the threshold for 
the structure functions is now defined in terms of v or T4^, rather than W. Hence, we consider two different cases in 
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choosing the £ integration limits for the ratio (|25| . First, for a given Q 2 , we choose the £ limits as in Eq. (|26|) . That 
amounts to defining them by the condition 

£min - £(VF = 1.6 GeV, Q 2 ), £ max = i{W = 1.1 GeV, Q 2 ) . (27) 

We refer to this choice as integrating "from 1.1 GeV". The integration limits for the DIS curve always correspond 
to this choice. As a second choice, for each Q 2 we integrate the resonance curve from the threshold, that is from as 
low W as achievable for the nucleus under consideration. This corresponds to the threshold value at higher £ and 
is referred to as integrating "from threshold" . With this choice we guarantee that the extended kinematical regions 
typical for resonance production from nuclei are taken into account. Since there is no natural threshold for the £ m i n , 
for both choices it is estimated from W — 1.6 GeV, as defined in Eq. (f2"T)) . 

The results for the ratio in Eq. (|2"5|) are shown in the right panel of Fig. [U The curve for the isoscalar free-nucleon 
case is the same as in Ref. [3| with the "GRV" parameterization for the DIS structure function. One can see that the 
carbon curve obtained by integrating "from threshold" lies above the one obtained by integrating "from 1.1 GeV", 
the difference increasing with Q 2 . This indicates that the threshold region becomes more and more significant, as one 
can see from Fig. [3J The closer the ratio (|2~5t gets to 1, the higher the accuracy of local duality is. Our calculations 
for a carbon target show that: 1) the ratio grows with Q 2 , just like in the isoscalar frcc-nuclcon case; 2) the ratio is 
lower than the free-nucleon value for both choices of the integration limits. This means that the integrated resonance 
contribution is always smaller than the integrated DIS one. In search for an explanation for this discrepancy, we 




FIG. 5: Electromagnetic structure functions F2 in DIS region for a free isoscalar nucleon as obtained via the GRV parameteri- 
zation at Q 2 = 10 GeV 2 (solid curve) and for a carbon nucleus as measured experimentally at Q 2 — 30, 45 and 50 GeV 2 . 

compare how nuclear effects influence the resonance and DIS curves. As it was illlustrated in Fig. [3J the nuclear 
effects suppress the resonance peaks by 40%-50%, broaden them and shift them to lower £ values. The experimental 
DIS values for the carbon nucleus, on the other hand, are only 5% — 10% lower than the DIS curve for the free isoscalar 
nucleon. This is illustrated in Fig. El where the DIS structure function F2 for a carbon target is compared to the 
GRV parameterization for the free isoscalar nucleon at Q 2 = 10 GeV 2 . In conclusion, one can say that nuclear effects 
have a much more dramatic effect in the resonance region than in the DIS regime. 
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Similar calculations can be done for other nuclei. First of all, it would be interesting to compare an isoscalar 
nucleus with a nucleus with neutron excess. We show the structure functions (F| 56Fe /56) and (F| 52fe /52) versus 
W for several Q 2 values in the left panel of Fig. [5] The structure functions for 52 Fe are only marginally higher than 
those for 56 Fe. This can be explained by the fact that the electromagnetic A-production cross section is equal for 
proton and neutron targets. In the second resonance region, the cross sections on the proton are typically 5% — 30% 
higher than those for the neutron. The overall effect, however, is hardly visible for an excess of 4 neutrons out of 56 
nucleons. From an experimental point of view, it is also interesting to compare carbon with iron target nuclei. In the 
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FIG. 6: (color online) (Left) Structure functions F2 for electroproduction on iron-52 and iron-56 versus W. Curves are for 
Q 2 = 0.2, 0.45, 0.85, 1.4, and 2.4 GeV 2 (indicated on the spectra). (Left) Ratio (Ff 6 Fe /5&) / (F% 12 c / 12) versus Bjorken variable 



x compared to DIS data [3ll 



right panel of Fig. [6j we plot the ratio of structure functions (F| 56 - Fe /56) / (F| c /12) versus x for several values of Q 2 
ranging from 0.2 GeV 2 to 3.3 GeV 2 . All curves are shown in the £ region corresponding to 1.1 GeV < W < 2.0 GeV. 

We stress that there is little physical meaning in the fine structure of the curves in the right panel of Fig. O The 
peaks in the curves, for example, do not coincide with the resonance peaks. As one can see, the iron structure 
functions appear to be very close to the carbon ones: for each Q 2 the ratio of the iron to carbon structure functions 
does not deviate more than 5% from the value of 1. When averaged, this ratio slightly decreases with increasing Q 2 , 
a behavior which is also exhibited by the DIS data presented in the same figure. The latter were measured by the 
NMC Collaboration 31j, the mean Q 2 in the experiment varying from 20 GeV 2 for x ~ 0.1 to 60 GeV 2 for x — > 1. 



IV. DUALITY IN NEUTRINOPRODUCTION 

In a previous paper [13I ] it was demonstrated that in neutrino reactions quark-hadron duality does not hold for 
proton and neutron targets separately. This is a principle feature of neutrino interactions, stemming from fundamental 
isospin arguments. For the charged current reaction v^p — ► ji" R ++ , only isospin-3/2 R ++ resonances are excited, in 
particular the P33(1232) resonance. Because of isospin symmetry constraints, the neutrino-proton structure functions 
for these resonances are three times larger than the neutrino-neutron ones. In neutrino-neutron scattering, both the 
isospin-3/2 resonances and the isospin-1/2 resonances contribute to the structure functions. The interplay between 
the resonances of different isospins allows for duality to hold with reasonable accuracy for the average over the proton 
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and neutron targets. It appears reasonable that one may expect a similar picture to emerge in neutrino reactions 
with nuclei. 
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FIG. 7: (color online) The computed resonance curves F% /56 and Xf^F^ /56 as a function of £, for Q 2 = 0.2, 0.45, 0.85 



at Q 2 DJ s = 7.94, 12.6 and 19.95 GeV 2 . 



1.4, and 2.4 GeV . The calculations are compared with the DIS data from Refs. [26, 271. The DIS data refer to measurements 
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FIG. 8: (color online) Ratios 7| Fe and if Fe defined in Eq. J25J for the free nucleon (dash-dotted line) and 56 Fe. For 56 Fe 
the results are displayed for two choices of the underlimit in the integral: W — 1.1 GeV (solid line) and threshold (dotted line). 
For each of these two choices we have used two sets of DIS data in determining the denominator of Eq. ([25} . These sets of DIS 
data are obtained at Q 2 DIS = 12.59 and 19.95 GeV 2 . 



The structure functions and xaF^ for neutrino-iron scattering are shown in Fig. [7] The curves for the isoscalar 



free nucleon case is identical to the one presented in Ref. 13| with the "fast" fall-off of the axial form factors for the 
isospin-1/2 resonances. 

Like for the electron-carbon results of Fig. [31 the resonance structure is hardly visible. Indeed for each Q 2 the 
computed resonance curves display one broad peak. The resonance structure functions are compared with the exper- 
imental data in DIS region obtained by the CCFR [26( and NuTeV [27j collaborations. It appears that the resonance 
curves slide along the DIS curve, which indicates global duality. Like for the electron results discussed in previous 
section, however, the resonance and xaF^ predictions are noticeably lower than the DIS measurements. 
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The ratios 1% Fe and if Fe denned in Eg. ([25]) are shown in Fig. [8] Our results show, that 1) these ratios are 
significantly smaller than 1; 2) they are significantly smaller than the one for the free nucleon ; 3) I2 is lower than 
the corresponding ratio for electroproduction; 4) I2 and I3 slightly decrease with Q 2 which is the opposite behavior 
of what was observed for electrons. 
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FIG. 9: (color online) (Left) Weak structure functions F2 in the resonance region for a free neutron (solid curve), for Is 1 / 2 
(dashed curve), lp 3 ^ 2 (dash-dotted curve), Id 5 / 2 (long-dashed curve) and lf 7 ^ 2 (short-dashed curve) neutrons in B6 Fe. The 
three sets of curves correspond to Q 2 = 0.2, 0.85, and 2.4 GeV 2 . (Right) Weak structure functions F2 in DIS region for a free 
isoscalar nucleon as obtained via GRV parameterization at Q 2 = 10 GeV 2 (solid curve) and for iron-56 nucleus as measured 
experimentally Q,H at Q 2 = 7 - 9 T 12.59 and 19.95 GeV 2 . 



In an attempt to explain the above observations, we compare the free isoscalar structure functions with the 56 Fe 
ones. In the left panel of Fig. [9] the structure function for a neutron in the Is 1 / 2 , lp 3 / 2 , Id 5 / 2 and l/ 7 / 2 iron 
shells are contrasted with the structure function for a free neutron. For a fixed Q 2 , the peak at the larger value of 
the Nachtmann variable corresponds to the A resonance. The peak at smaller £ corresponds to the second resonance 
region. It is clear that the nuclear effects reduce the peaks by about 30 — 50% and shift them to lower £ values in 
comparison with the free nucleon case. The suppression is most significant for the single-particle shells close to the 
Fermi surface. 

In close resemblance to what was observed in the discussion of the electron-nucleus cross sections of previous section, 
the peculiar Fermi smearing pattern for each shell introduces additional averaging when summing over shells. For 
a bound proton, the effect of suppression is nearly the same. In the right panel of Fig. [9l measured DIS structure 
functions for 56 Fe at various Q 2 are compared to the GRV parameterization (Q 2 = 10 GeV 2 ) for a free isoscalar 



nucleon 



13l | . It is obvious that the measured nuclear DIS structure functions are very similar to the free- nucleon ones. 



Thus, we predict a substantial nuclear reduction of the resonance strength, whereas the data in the DIS region do not 
point to such a reduction. This explains the computed low values of the ratios in Fig. [HJ 

We wish to stress that the low values of lf BFe and lf &Fe are not related to the neutron excess. We remind that 
the neutron structure functions for the A resonance are 3 times smaller than the proton ones. The structure function 
for isoscalar 52 Fe is only about 5% larger than for 56 Fe. This is shown in Fig. [10] for F2 and xaF^. The effect can be 
easily estimated from (26 • 3 • / + 30 • /)/ (26 • 3 • / + 26 • /) « 1.04, where / is the neutron structure function in the A 
region. In the second resonance region the difference is even smaller. 
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FIG. 10: (color online) Structure functions (left) and xaF^/A (right) for neutrinoproduction on iron-52 and iron-56 versus 
W. Curves in the resonance region are for Q 2 = 0.2, 0.45, 0.85, 1.4, and 2.4 GeV 2 (indicated on the spectra). 
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FIG. 11: (color online) Ratios {F%™ Fe /m)/{F% 12a /12) (left) and {xse Fe F^ 56pe/56) / (x 5 2 Fe F% 12 ° /12) (right) versus Bjorken 
variable x for Q 2 = 0.2, 0.45, 0.85, 1.4, 2.4 and 3.3 GeV 2 



It is also interesting to make a comparison with the carbon nucleus. The ratios of iron to carbon structure functions 
F% and xaF^ versus x are shown in Fig. QT] For each Q 2 the £ range corresponds to 1.1 < W < 2.0 GeV. Like in 
the case of electromagnetic reaction, the ratios are close to 1, but a bit lower in general and the average is slightly 
increasing with Q 2 . Remark that the peaks in Fig. [TT] are not related to resonances and that the fluctuations which 
are of the order of 5% can be attributed to subtleties in the shell structure of the various target nuclei. 



V. SUMMARY 

In view of the current experimental activities, there is great need for an efficient framework for reliably predicting 
neutrino nucleus cross sections and for a deeper understanding of quark-hadron duality in nuclei. We performed a 
phenomenological study of duality in electron-nucleus and neutrino-nucleus structure functions. 

Using the Dortmund-group model for the production of the first four lowest-lying nucleon resonances and using 
single-particle wavefunctions from the Hartree approximation to the relativistic aw model, we computed the structure 
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functions xaF^, F^ and xaF^ in the resonance region for carbon and iron targets and compared them with the 
measured DIS ones. At the same time we compared the computed resonance structure functions for nuclei with those 
for a free nucleon. For quantitative comparisons, we defined the ratios h{Q 2 ) of integrated resonance to DIS structure 
functions. Perfect quark-hadron duality is reached for I%{Q 2 ) values of unity. 

Summarizing our results, we observe that the computed resonance contribution to the lepton-nucleus structure 
functions is qualitatively consistent with the measured DIS structure functions. This means that global quark-hadron 
duality holds for nuclei. The computed integrated resonance strength, however, is about half of the measured DIS 
one. Contrary to the free nucleon case, where the ratios h{Q 2 ) are at the level of 0.8, we find for nuclei 0.6 for 
electroproduction and 0.4 for neutrinoproduction. This points towards a scale dependence in the role of the nuclear 
effects. It is obvious that nuclear effects act differently at lower Q 2 (resonance regime) than at higher Q 2 (DIS regime). 

In our presented analysis we include the resonance contributions and ignored the role of the background terms. 
Further investigations require a theoretical or phenomenological model for the background contributions in the first 
and second resonance region. One could for example estimate the role of the background contribution to the A- 
resonance region within the context of the non- linear sigma model [32| . Extending these or similar models to higher 
W values and incorporating them in a model for lepton reactions with nuclei could be the next step in exploring 
quark-hadron duality. 
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